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1  Statement  of  the  problem  studied  and  Sig¬ 
nificance 

This  proposal  concerns  development,  analysis,  implementation,  and  appli¬ 
cations  of  efficient  and  accurate  numerical  methods  for  interface  problems, 
especially  with  moving  interface  and/or  free  boundary  problems,  and  prob¬ 
lems  defined  on  irregular  and  infinite  domains.  Specifically,  we  will  develop 
numerical  methods  for  fully  elastic  systems  with  composites  and  material  in¬ 
terfaces,  fully  transient  simulations  of  two  phase  flows,  e.g.,  droplets  deforma¬ 
tion  and  breakups  in  shear  flows  of  immiscible  fluids,  and  electro-magnetics, 
elastic  wave  propagation  in  heterogeneous  media. 


2  Summary  of  the  most  important  results 

•  We  have  developed  a  fast  iterative  solver  for  scattering  by  elastic  objects 
in  layered  media. 

•  We  have  nice  convergence  analysis  of  the  primal-dual  active  set  strategy 
for  diagonally  dominant  systems. 

•  We  have  developed  a  preconditioned  iterative  methods  on  sparse  sub¬ 
spaces. 

•  We  have  developed  new  immersed  finite-element  methods  for  elliptic 
and  elastic  interface  problems  with  homogeneous  and  noil-homogeneous 
jump  conditions  [2,  3].  The  new  methods  do  not  require  body-fitted 
meshes.  Non-homogeneous  jump  conditions  are  removed  by  a  known 
function  that  has  the  same  jump  conditions  as  that  of  the  solution  via  a 
level  set  function.  The  new  approaches  provide  a  second  order  discrete 
delta  function  for  elliptic  and  elastic  interface  problems.  The  results 
have  been  published  in  SIAM  Numer.  Anal.  [3]. 

•  In  the  past,  we  derived  jump  conditions  for  the  Stokes  and  Navier- 
Stokes  equations  with  continuous  viscosity  [7,  8].  We  have  derived  the 
jump  conditions  for  discontinuous  viscosity  recently  using  the  idea  of 
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the  immersed  interface  method  in  [11],  We  also  developed  the  immersed 
interface  methods  for  various  problems  including  new  approach  to  deal 
with  pressure  jump  conditions  [13],  the  Stokes  flows  with  discontinu¬ 
ous  viscosity  [10,  18],  biharmonic  equations  on  irregular  domains  [1], 
incompressible  Navier-Stokes  equations  on  fixed  irregular  domains  [10] 
with  application  to  flow  past  fixed  obstacles. 

•  In  the  application  to  problems  in  mathematical  biology,  our  immersed- 
interface/level  set  method  applied  to  the  biological  problem  of  forces 
creating  branching  morphogenesis  shows  that  contractility  of  the  mes¬ 
enchyme  is  indeed  sufficient  to  create  a  realistic  cleft  in  a  branching 
embryonic  rudiment  of  a  glandular  organ.  This  is  quite  possibly  the 
first  simulation  ever  to  realistically  model  the  forces  in  branching  mor¬ 
phogenesis,  [15]. 

•  We  have  some  detailed  analysis  for  a  mixed  elliptic  and  hyperbolic 
system  for  modeling  tissue  deformation  [6]. 

With  partial  support  of  the  ARO  grant,  we  have  published  about  15 
journal  papers  [5,  13,  10,  3,  11,  18,  1,  15,  14,  6,  4,  20,  12,  16,  19,  17], 
and  one  book  [9]  in  the  last  three  years. 

3  Other  contributions 

Please  see  the  list  of  publications  of  and  other  reported  information  in  the 
on-line  form  of  the  final  report. 

In  a  summary,  we  are  able  to  accomplish  almost  all  the  projective  in  our 
original  proposal  and  more.  Furthermore,  the  ARO  grant  really  helped  us  to 
get  our  research  projects  going  and  move  to  new  research  areas.  We  really 
appreciate  ARO  for  the  support. 
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